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ABSTRACT 

In this paper we present a detailed spectroscopic analysis of the suspected marginal 
Am star HD 71297. Our goal is to test the accuracy of two different approaches to 
determine the atmospheric parameters effective temperature, gravity, projected rota- 
tional velocity, and chemical abundances. The methods used in this paper are: classical 
spectral synthesis and the Versatile Wavelength Analysis (VWA) software. 

Since our star is bright and very close to the Sun, we were able to determine its 
effective temperature and gravity directly through photometric, interfcrometric, and 
parallax measurements. The values found were taken as reference to which we compare 
the values derived by spectroscopic methods. Our analysis leads us to conclude that 
the spectroscopic methods considered in this study to derive fundamental parameters 
give consistent results, if we consider all the sources of experimental errors, that have 
been discussed in the text. In addition, our study shows that the spectroscopic results 
are quite as accurate as those derived from direct measurements. 

As for the specific object analyzed here, according to our analysis, HD 71297 
has chemical abundances not compatible with the previous spectral classification. We 
found moderate underabundances of carbon, sodium, magnesium, and iron-peak ele- 
ments, while oxygen, aluminum, silicon, sulfur, and heavy elements (Z 39) are solar 
in content. This chemical pattern has been confirmed by the calculations performed 
with both methods. 
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1 INTRODUCTION 

The current era is characterized by an enormous growth of 
stellar data of different nature. The advent of space mis- 
sions, suc h as CoRoT ( Conve ction, Rotation and planetary 
Tran sits; iBaglin et all I 20071 ) and Kepler (|Borucki et al.l 
1 19971 ). designed to obtain precise photometry for an impres- 
sive number of stars, have led astronomers to undertake sev- 
eral ground-based spectroscopic campaigns in order to get as 
accurate as possible fundamental parameters such as effec- 
tive temperature, surface gravity, projected rotational veloc- 
ity, and metallicity. The analysis of this huge amount of data 



* Based on observations made with the Italian Telescopio 
Nazionalc Galileo (TNG) operated on the island of La Palma by 
the Fundacin Galileo Galilei of the INAF (Istituto Nazionale di 
Astrofisica) at the Spanish Obscrvatorio del Roque de los Mucha- 
chos of the Instituto de Astrofisica de Canarias 
f E-mail: gca@oact.inaf.it 



is possible, on reasonably short timescales, only through the 
use of automatic or semi-automatic procedures. 

Among the various targets observed by Kepler, Am 
stars are assuming an ever growing importance. In fact, it 
was once thought that Am stars did not pulsate, and that 
the explanation for this is that atomic diffusion is expected 
to drain helium from the Hen driving zone. More intensive 
observations have revealed that this picture is not correct 
and several Am stars are known to pulsate from ground- 
based observations (based on superWASP photometry, see 
Smallev et afl|201ll ) as well as from Kepler satellite data (see 



Balona et al.ll201ll ). In this last paper, the authors showed 



that an accurate measure of the location of the pulsating 
Am stars in the HR diagram is crucial to pin down effec- 
tively the pulsation models and to put observational con- 
straint on the instability strip for pulsation in Am stars. 
However, for most of the investigated objects, a modern de- 
termination of the stellar parameters such as effective tem- 
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peratures, gravities and chemical abundances, based on high 
resolution spectroscopy is still lacking. To fill this gap, we 
have already started a spectroscopic campaign devoted to 
obtain high resolution spectra of Am stars, being the first 
data obtained with the spectrograph SARG installed at the 
Italian telescope TNG. Our purpose in the near future is 
to analyze this growing amount of spectra to derive accu- 
rate stellar parameters for the investigated objects in the 
shortest possible time. 

There are several examples in the literature of 
papers regarding various methods for aut omatic (or semi- 
autom a tic) spectrosco pic ana lysis, ( Tkachcnko et all 



20121; iLehmann et all l201ll ; iNiemczura et al.1 |2009| : 
Erspamer fc North! 12003 ). However, it is important to 
ascertain the accuracy of the derived parameters obtained 
using different methods. A few recent papers have partially 
addressed this problem, reporting that significant differences 
can be found both in the derived stellar parameters (see e.g. 
iFossati et al.l l201ll ) and/or in the estimated uncertaintie s 
in these quantities (see e.g. iMolenda-Zakowicz et alj|201ol ). 

With the aim of understanding how accurately we can 
derive the stellar parameters and the HR location for our 
sample of Am stars, we decided to concentrate our efforts 
on comparing the results that can be obtained by analyz- 
ing the high resolution spectrum for one star in our sample 
with the two methods available to u s : i) classical spectral 
synth esis (see ICatanzaro et al.l l201ll ; ICatanzaro fc Balonal 
12012] , and references therein) and ii) Ver satile Wavelength 
Analysis (VWA) fsee lBruntt et al.1 l2010aH bl. and references 
therein) . 

To carry out this exercise, we choose to analyze in detail 
the Am star HD 71297. This is the most suitable object in 
our sample because it is very bright (V = 5.58 mag) and 
the SARGQTNG spectrum has a very good signal-to-noise 
ratio. Moreover, HD 71297's brightness allowed us to find 
very useful observational data in the literature, including an 
accurate Hipparcos parallax and interferometric measures 
which allow us to estimate independently from spectroscopy 
the stellar parameters. Thus, the main goal of this paper is to 
ascertain if the two quoted methods of quantitative spectral 
analysis lead to different or, on the contrary, to comparable 
results when used to exploit exactly the same observational 
material. 



This star has also been studied bv lGuthriei 1 19871 ) who 
derived atmospheric parameters, T c fi =7900 K, logg=4.2, 
and calcium abundance, log iV(Ca) = 6.33, expressed in a 
scale in which l og N(H) = 12.0, that is p ractically the solar 
value. Later on, iKuenzli fc North! (|l998l ) revisited the star, 
finding it a little bit cooler and less evolved, T c ff =7712 K 
and log g=4.06, but still with solar calcium abundance. 

No other studies have been found in the recent litera- 
ture, especially no detailed analysis of the chemical pattern 
of its atmosphere. Thus, we chose this object as a benchmark 
for the two methods of analysis we would like to compare in 
this paper. 



3 OBSERVATION AND DATA REDUCTION 

Spectroscopic observations of HD 71297 were carried out 
with the SARG spectrograph, which is installed at the Tele- 
scopio Nazionale GaMeo, located in La Palma (Canarias 
Islands, Spain). SAR G is a high-resolutio n cross-dispersed 
echelle spectrograph (jGratton et al.|[200ll ) that operates in 
both single-object and longslit (up to 26) observing modes 
and covers a spectral wavelength range from 370 nm up to 
about 1000 nm, with a resolution ranging from R = 29 000 
to 164 000. 

Our spectra were obtained on 2011, February 21 at 
R = 57 000 using two grisms (blue and yellow) and two fil- 
ters (blue and yellow). These were used in order to obtain a 
continuous spectrum from 3600 A to 7900 A with significant 
overlap in the wavelength range between 4620 A and 5140 
A. We acquired the spectra for the star with an exposure 
time of 120 sec. and a signal-to- noise ratio S/N of at least 
100 per pixel in the continuum. For example, in the region 
centered around the Mgl triplet at 5170 A the measured S/N 
is about 120. 

The reduction of all spectra, which included the sub- 
traction of the bias frame, trimming, correcting for the 
flat-field and the scattered light, the extraction for the or- 
ders, and the wavelength calibration, was done using the 
NOAO/IRAF package^]. The IRAF package rvcorrect was 
used to make the velocity corrections due to Earth's motion 
to transform the spectra to the bary centric rest frame. 



2 THE TARGET. THE AM: STAR HD 71297 

H D 71297 has bee n classified as suspected marginal Am star 
by ICowlevI (|l968l ). with a spectral type of A5m: (colon de- 
noted its marginal nature) . We remind readers the definition 
of marginal Am star as an Am star in which there is a dif- 
ference of less than five subclasses between the K-line and 
the metal line spect ral types , and i n which the line strength 
anomalies are mild. ICowlevI (|l968h concluded that the star 
exhibits only weak metallic lines, looking like the spectrum 
of P Ar i. This classification has been confirmed one year 
late r bvlCowlev et al.l (Il969l ) . 

lAbt fc Levvl (| 19851 ) searched for binar ty among a sam- 
ple of 60 Am stars and concluded that HD 71297 may 
be variable in radial velocity with a period of a hundred 
days. They assigned the spectral type of kA8hA9mF0 to 
HD 71297. 



4 STELLAR PARAMETERS FROM 

PHOTOMETRY, INTERFEROMETRY AND 
PARALLAX. 

Before starting with the analysis of the spectrum of 
HD 71297, it is extremely important to estimate the stel- 
lar parameters, mainly T e fi and logp, through independent 
methods, such as those based on intermediate-band photom- 
etry, parallax, interferometry etc. This is useful to have: i) a 
reliable starting point for the spectroscopic analysis and ii) 
a trustworthy and accurate independent evaluation of the 
stellar fundamental parameters to compare with. 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc. 
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Table 1. Evaluations of T e jj and log g values of HD 71297 on the basis of various photometric systems (see text for details). 



Photometry 



L off 



log 5 



Calibration 



uvby-fi 
uvby-f3 

UBVB 1 B 2 V 1 G 

uvby-fi 
uvby-fi 
uvby—f3 



7700±40 
7970±30 
7770±80 
7770±125 
7800±125 
7780±125 



4.06±0.04 
4.28±0.04 
4.46±0.06 
4.08±0.125 
4.06±0.125 
4.06±0.125 



Napiwotzki et al. (1993) 

Ribas et al. (1997) 
Kuenzli fc North (1998) 
Smallev fc Kupka (19971 
Heiter et al. (2002') MLT o=0.5 



Hcitc r et al. (20021 using Canuto. Goldman fc M azzitclli (1996) 



4.1 Evaluation of T e g and log g from Stromgren 
and Geneva photometry 

A first estimate of T G ff and log g for HD 71297 can be ob- 
tained from the Stromgren photometry: V = 5.607 ± 0.002, 
b~y = 0.123 ±0.002, mi = 0.197±0.004, Cl = 833 ±0.006, 
/3 = 2.831 ± 0.003 (|Hauck fc Mermillodl fl99Sl ). However, 
It is important to check the value of V^-band magnitude 
because HD 71297 is a (low-amplitude) variable star (see 
iBalona et ail 120111 ). A safe average V magnitude can be 
derived from the time-series photometry in the Hippar- 
cos (Hp) and Tycho (Vr) photometric systems. Includ- 
ing the uncer tainties in th e trasformations to the John- 
son V system (|BesseUll20o"ol ), we get V=5.59±0.01 mag and 
V r =5.60±0.01 mag. Finally, an additional ind ependent mea- 
sure was reported bv lKuenzli fc North! (|l998l ) who measured 
V=5. 604±0. 008 mag. On this basis, we adopted a value of 
V r =5.60±0.01 mag for HD 71297. This value is consistent 
within the errors with all the above quoted evaluations. 

The other Stromgren indices should be less uncertain 
with respect to V-band, however it is likely that the errors 
are somewhat underestimated. 

We adopted an u pdated version of the TempLogG^ 
software (1 Rogers! Il995l ) to estimate T c ff and log g by us- 
ing the calibrations present in the pac k age. W e d id not 
consid e r the older calibration by | Balonal (|l984j ) and iMoonl 
1 19851 ); iMoon fc Dworetskvl dl985l), b u t only the more re - 
cent ones~by iNapiwotzki et all (|l993T ); iRibas et all (|l997l ); 
iKuenzli fc North! ( 19981 ). Note that the last calibration is 



for the Geneva UBVB1B2V 1G system, wh ose photometry 
for HD 71297 is provided by iRufenerl (|l988l ). These results 
are summarized in the first three rows of Table U In addi - 
tion , we consider e d the results bv lSmallev fc Kupka! (|l997l ) 
and iHeiter et al.l (|2002i ) who provided uvby grids based on 
the Kurucz model atmospheres but with different treatment 
of the convection. In particular we report on Table[pthe T e g 
and lo gff estimated by interpolating the Smalle v fc Kupka! 
1997 ) grids that were built using the lCanuto fc Mazzitellil 



19911 ) convection treatment. Similarly, the last two col- 



umn of Table [T] list the stellar p arameters ob t ained us- 
ing two choices for the gridfl by IHeiter et al.l (|2002l ): i) 
standard mixing- lenght theory (MLTjf with a=0.5; ii) the 
ICanuto. Goldman fc Mazzitellil (|l996h treatment of the con- 
vection. In all these cases the error associated to the parame- 



ters was imposed as half-the width of the grid. An inspection 
of Table \T\ reveals some dispersion among the various cali- 
brations both in T ff and loggr, even if the last three values 
are absolutely equivalent within the errors. Quantitatively, 
a simple average of the results gives: (T c h)=7800±90 K and 
(log5)=4.17±0.17 dex. 

As for the interstellar reddening, TempLogG gives as a 
result E(b - y)=0.009±0.003 mag. This value is probably 
non-significant, because of the very low value of both red- 
dening and associated uncertainty (underestimated in our 
opinion). Since the star is only 50 pc far from the Sun (see 
next sub-section), it is acceptable to neglect the effect of the 
interstellar absorption. This assumption was also confirmed 
a posteriori by looking into the spectrum of HD 71297 for the 
interstellar Nal lines at 5890.0, 5895.9 A which can be used 
to est imate the interstellar reddening (see lMunari fc ZwitteJ 
Il997l ). Indeed, these lines are practically undetectable. To 
conclude, in the following analysis we have neglected the 
interstellar reddening. 

4.2 Evaluation of the stellar parameters from 
Photometry, Parallax and Interferometry 

The Hipparcos para llax for HD 71297, n = 19.99 ± 0.38 mas 
( van Leeuwenll2007l ). allows us to estimate with the excellent 
accuracy of 2% the distance of this object: d = 50± 1.0 pc. 
It is then straightforward to calculate the visual absolute 
magnitude of the target: M v = 2.105 ± 0.045 mag. To de- 
rive the luminosity we need to evaluate first the visual bolo- 
me tric correction BCy. To this a im we adopted the models 
bv lBessell. Castelli fc Pled (( 19981 ) where it is assumed that 
M b oi,0 = 4 74 ma£|. We then get BC V = 0.03 ± 0.01 mag 
for HD 7129 7. We note that to obt a in thi s value we inter- 
polated the iBessell. Castelli fc Pled (|l998l ) tables by using 
as input the T e ff and logg obtained from the photometry. 
The error was conservatively assumed to be 0.01 mag to 
take into account the observed dispersion of the photomet- 
rically derived T e fi and log g. We have now estimated all the 
quantities needed to calculate the luminosity of HD 71297: 
log(L/L ) = 1.04±0.0fl 

Another very important piece of information is repre- 
sented by the angular diameter of HD 71297, which was pub- 
lished by lLafrasse et all (|2010l ): a = 0.367 ± 0.025 mas. By 
using this value, together with the distance derived from 



2 available through http://www.univie.ac.at/asap/manuals/tipstricks/tcmlj^^ uncertainty on the bolometric magni- 

3 These grids are available on the NEMO site tude of the Sun, in all the calculations we associated an error of 
www.univie.ac.at/nemo/gci-bin/dive.cgi 0.01 mag to this value 

4 defined as the ratio a=l/H p of convective scale length / and 6 Whenever possible, we have summed in quadrature the uncer- 
local pressure scale height H v tainties on the single quantities when propagating the errors 
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Table 2. Astrophysical quantities for HD 71297 evaluated inde- 
pendently from spectroscopy. 



Parameter 



Value 



a 

My (mag) 
M b ol (mag) 
T off (K) 

log(L/L ) (dex) 
log g (dex) 
R/R 
M/M Q 
age (Myrs) 



19.99 ±0.38 
0.367±0.025 
2.105 ±0.045 
2.14 ±0.04 
7810 ± 90 
1.04 ±0.02 
4.17 ±0.06 
1.97 ±0.14 

i 77 +0.12 
1,1 '-0.19 
790 ± 90 




7500. 7000. 
Effective temperature (K) 



Hipparcos parallaxes, d = 50 ± 1 pc, it is straightforward to 
evaluate the radius of our target: R= 1.97±0.14 Rq. 

With the observational data available for HD 71297, the 
most convenient way to estimate the effective temperature 
is from the definiti on of surface brightness as given in Eq. 
14.8 of lGravl <|2005h . After simple algebra, we get: 

log Toff = c- 0.1 • V + 0.1 • BC- 0.5 log (1) 

where On is the angular radius in arcsec and c is a constant 
given by: 

c = log T e | + 0.1 • m v - 0.1 ■ BC Q + 0.5 log 9% (2) 

To est imate c we used the Sun "constants" re- 
ported by iBessell. Castelli fe Pled (|l99Sf ): T o ® =5781±4 
K, m® = -26.76+0.01 mag, BCq = -0.07+0.01 mag, and 
0% = 959.61+0.01 arcsec. Thus c = 2.58405T0.00140, and, 
from Eq. [TJ T eff = 7850 ± 280 K. In the error budget the 
most important contributor is the angular radius (diame- 
ter) at the ~3.5% level, being the contribution of the other 
quantities practically negligible. The effective temperature 
value derived here is almost indistinguishable from that es- 
timated directly from Stromgren and Geneva photometry, 
but with a larger uncertainty. Nevertheless, this evaluation 
appears robust, since it is based on different kind of inde- 
pendent measures (photometry, interferometry). Moreover, 
also the uncertainty can be calculated in a reliable and phys- 
ical way. Hence, we decided to take a weighted average of 
the effective temperatures obtained with the two methods 
as the value to be used as a reference for comparison with 
spectroscopy. As a results, we obtain T a g = 7810 ± 90K. 

To estimate logp, it is convenient to use the following 
expression that can be easily derived from the definition of 
g and from the Stefan-Boltzmann relation: 

logg = 4 lo g (T 8fl /T%) + log(M/M ) + 2 log(Tr) 

+0.4(1/ + BC V + 0.26) + 4.44 (3) 

Where the various terms of the equation have the usual 
meaning and M/M is the ratio mass of the star over mass 
of the Sun. Before using Eq. [3] we have to evaluate the mass 
of HD 71297. This can be safely done by adopting the cal- 
ibration mass -My (valid for luminosity class V stars) by 
iMalkovl (|2003l ) that was derived on the basis of a large sam- 
ple of eclipsing binaries stars. Hence, by using our My es- 
timate, we evaluated log(M/M ) = 0.248+ 0.05 dex (or 



Figure 1. Location of HD 71297 in the HR diagram together with 
evolutionary tracks and isochrones for logf ranging from 8.85 to 
8.95 (step 0.05 and t in yrs). 

M/M® = 1.77±°;^), where the error is largely conservative 
and completely dominated by the dispersion of the mass- 
My relation. 

Finally, applying Eq. [3] we obtain log g — 4.17 ±0.06, 
again, in excellent agreement with the purely photometric 
evaluation but with a much smaller uncertainty. We will 
then use this value for comparison with spectroscopy. 

Before closing this section, for completeness, it is worth 
estimating the age of HD 71297, taking advantage of the 
logTcff and log(L/L0) values estimated in this section. 
The location of the star in the HR diagram, together with 
some evolutio nary tracks comput ed for the solar metallic- 
ity Z = 0.019 (|Girardi et alJl200Ch . is shown in F i g. [T] Also 
shown are isochrones computed bv lMarigo et alj (|2008l ) for 
the same Z (=0.019) and for ages ranging from log t = 8.85 
to 8.95, in steps of 0.05 (t in years). The location of the 
star indicates a mass of M ~1.75 ± 0.05 Mq (in perfect 
agreement with the value adopted above) and an age of 
t ~ 790 ± 90 Myrs. 

In Table [2] we summarized all the astrophysical quan- 
tities for HD 71297 evaluated independently from spec- 
troscopy. 



5 ATMOSPHERIC PARAMETERS FROM 
SPECTROSCOPY 

In this section we will present separately the results from 
the analysis of the high resolution spectra of HD 71297 as 
obtained with the two different approaches quoted in the 
introduction. We want to stress that the two analysis have 
been performed by using the same spectrum, prepared as 
described in Section [3] In this way, we can be sure that any 
differences arising during the analysis will not be attributed 
to the quality of the data, but directly to the method itself. 

5.1 Stellar parameters and abundance analysis 
with classical spectral synthesis 

This a pproach was already discus sed in ICatanzaro et al] 
l|201ll) : ICatanzaro fc Balonal (|2012l ) and references therein. 
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Figure 2. Observed Hp and with over-imposed the synthetic 
ones computed with the fundamental parameters derived in this 
section. 

Thus, here we briefly summarize the principal features of 
this method. In order to determine the optimal parameters, 
we minimize the difference between the observed and syn- 
thetic spectrum. Thus we minimize 

2 _ _j_ V - ^ / Jobs — ^th \ 2 
x - N 2^ V 5/obs J 

where N is the total number of points, I ODB and I t h are 
the intensities of the observed and computed profiles, re- 
spectively, and Slobs is the photon noise. Synthetic spectra 
were generated in three steps. First, we com puted a LTE at - 
mospheric model using the ATLAS9 code |Kurucdll993ah . 
The stellar spectrum wa s then synthesized using SYNTHE 
(|Kurucz fc Avrettl Il98ll ). Finally, the spectrum was con- 
volved with the instrumental and rotational profiles. 

As starting values of T e ff and log g, we used the values 
derived in the previous section. 

To decrease the number of parameters, we computed the 
vsini of HD 71297 by matching synthetic line profiles from 
SYNTHE to a number of metallic lines. The Mgl triplet 
at AA5167-5183 A was particularly useful for this purpose. 
The best fit was obtained for v sin i = 7.0 ±0.5 kms" 1 . This 
value is slightly lower than the one published bv lRover et al.l 
|2002l ) of usini = ll km s _1 . 

Uncertainties in T e g , log g, and v sin i were estimated by 
the change in p arameter values whic h leads to an increases 
of x 2 by unity (|Lampton et al.ll 19761 ). 

To determine stellar parameters as consistently as pos- 
sible with the actual structure of the atmosphere, we per- 
formed the abundance analysis by the following iterative 
procedure: 

(i) T e ff is estimated by computing the ATLAS9 model 
atmosphere which gives the best match between the ob- 
served H/3 and lines profile and those computed with 




5170 5175 5180 5185 5190 

Wavelength (A) 



Figure 3. Observed Mgl triplet with over-imposed the synthetic 
one computed with the fundamental parameters derived in this 
section 

SYNTHE (see Fig. [2j . The model was computed using solar 
opacity distribution function (ODF) and microturbulence 
velocity £ = 2.4 ±0.5 km s _1 , the latter value has been cal- 
cu lated following the cal i bratio n g = £(T ff , log g) published 
bv lAllende Prieto et al.l (|2004h . These two lines are the only 
useful for our purpose since they are located far from the 
echelle orders edges so that it was possible to safely recover 
the whole profiles. The simultaneous fitting of two lines led 
to a final solution as the intersection of the two \ 2 i so ~ 
surfaces. Another source of uncertainties is due to the diffi- 
culties in normalization as is always challenging for Balmer 
lines in echelle spectra. We quantified the error introduced 
by the normalization in at least 100 K, that we summed in 
quadrature with the error obtained by the fitting procedure. 
The final result is: 

T cff = 7500 ± 180 K 

The surface gravity was estimated from line profile fitting of 
Mgl lines with developed wings. This method is based on the 
fact that the wings of the Mgl triplet at AA 5167, 5172, and 
5183 A lines are very sensitive to log g variations. In practice, 
we have first derived the magnesium abundance through the 
narrow Mgl lines at AA 4571, 5528, 5711 A and the Mgn 
A 7877 A, and then we fitted the wings of the triplet lines 
by fine tuning the logg value. The magnesium abundance 
we got from these lines was log Alg/N to t = —4.62 ± 0.14. 
To derive \ogg by fitting spectral wings is essential to take 
into account very accurate measurements of the atomic 
parameters of the transitions, i.e. log gf and the radia- 
tive, Stark and Van der Waals d amping constan t s. Re - 
garding log gf we used the values of lAldenius et al.l l|l997h . 
Van der Waals damping constant is that calculated by 
iBarklem. Piskunov fc O'Maral (|200Ch (log 7Waals = —7.37) , 
the Stark damping constant is from iFossati et al.l (|201 lh 
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Figure 4. The figure shows six Fen lines in HD 71297 fitted by VWA (continuous line). The wavelengths of the fitted lines are given in 
the bottom right corner of each panel. 



(log7stark = —5.44), and the radiative damping constant 
is from NIST database (log7 rad = 7.99). 
This procedure results in the final value of: 

logp = 4.0 ± 0.1. 

(ii) As a second step we determine the stellar abundances 
by computing a synthetic spectrum that reproduced the ob- 
served one. Therefore, we divide our spectrum into several 
intervals, 25 A wide each, and derived the abundances in 
each interval by performing a x 2 minimization of the dif- 
ference between the observed and synthetic spectrum. The 
minimization algorithm has been written in IDL language, 
using the amoeba routine. W e adopted lists of spectra l lines 
and atomic parameters from lCastelli fc Hubrid (|2004|). who 
updat ed the parameters listed originally by Kurucz fc Belli 
(|l995h . 

For each element, we calculated the uncertainty in the 
abundance to be the standard deviation of the mean ob- 
tained from individual determinations in each interval of 
the analyzed spectrum. For elements whose lines occurred 
in one or two intervals only, the error in the abundance was 
evaluated by varying the effective temperature and gravity 
within their uncertainties given in Table [3] [T E g ± 5T e g] and 
[log g ± S log g] , and computing the abundance for T e g and 
logg values in these ranges. We found a variation of ~ 0.1 
dex due to temperature variation. We did not find any signif- 
icant abundance change by varying log g. The uncertainty in 
the temperature is the main error source on our abundances. 

T he derived abundanc es, expressed in terms of solar 
values jGrevesse et aLlfeoid ). are shown in Fig. [6] (red dots). 
For all the chemical elements for which we detected spectral 
lines in our spectrum, we found moderate under-abundances 
with respect to the solar analogues, with the exception of 
the elements with high Z (^ 40), for which we found normal 
abundances. 

As an example of our spectral synthesis we show in 
Fig. [3] the comparison between the observed spectral range 
from 5160 A to 5190 A and the synthetic one. In the second 
column of Table O we listed the abundances for our star 
derived with the method described above. 



5.2 Stellar parameters and abundance analysis 
with VWA 

The software package VWA relies on the iterative fit of syn- 
thetic profile regions around reasonably isolated absorption 
lines. VWA has a graphical user interface (GUI), which al- 
lows the user to investigate the spectra in detail, pick lines 
manually, inspect the quality of fitted lines, etc. 

VWA uses data from various sources. In particular, for 
th e calculation of synth e tic sp ectra it uses the SYNTH code 
bv lValenti fc Piskunovl ljl996T) . which works with ATLAS9 
models, and atomic parameters a nd line-broadening coef- 
ficients from the VALD database |Kupka et all 1 19991 ). We 
note that we used model atmospheres interpolated in the fine 
grid published bv lHeiter et aL d2002l) . These models rely on 
the original ATLAS9 code by Kurucz dl993a ]but use a more 
advanced convection description (iKupkal 19961 ) based on the 
works bv lCanuto fc Mazzitelli (|l992l ). 

A thorough de script io n of how VWA works can b e foun d 
in iBruntt et all (|2004l ); iBruntt. De Cat fc Aertsl (|2008h . 
Here we only recall the main steps of the analysis. 

• Normalization: the first step consists in an accurate nor- 
malization of the spectrum. This can be achieved by adopt- 
ing a synthetic spectrum (with approximately the same stel- 
lar parameters as the target) as a template to properly iden- 
tify (in the GUI) suitable regions to anchor the continuum 
of the observed spectrum. 

• Setup and line selection: atomic parameters and a pre- 
liminary model are setup on the basis of initial values for the 
stellar parameters. The line selection can be carried out ei- 
ther in automatic or manual way. We first searched for good 
lines (i.e. with low blending) in an automatic way, then we 
checked the result visually line by line on the GUI. 

• Fit of the lines and check of the result: each line is fitted 
by iteratively changing the abundance to match the equiv- 
alent width (EW) of the observed and calculated spectrum. 
The fitted lines are inspected in the GUI, problems with 
the continuum level or asymmetries in the line are readily 
identified, and these lines are discarded. This is done auto- 
matically by calculating the \ 2 of the fit in the core and the 
wings of the lines. This is followed by a manual inspection 
of the fitted lines. Figure [4] displays a few selected lines and 
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Figure 5. Abundance of Fe in HD 71297 for two different at- 
mosphere models as a function of equivalent width (left panels) 
and excitation potential (right panels). Fel and Fell lines are plot- 
ted with open and solid points, respectively. The horizontal line 
in each panel show the value of log NF e /N to ta i=-4.53 adopted 
for the Sun, according to lGrevesse et ah! l|2010f ). Top panels show 
the results for the adopted (best) parameters, while bottom pan- 
els display the effect of decreasing logg by 0.4 dex (see labels). 
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Figure 6. Chemical pattern for HD 71297 derived with the two 
different methods, red dots are from spectral synthesis and black 
squares from VWA. The bottom panel shows the difference be- 
tween the two methods, for all the chemical elements in common. 



the relative fits to demonstrate the quality of the observed 
data (narrow lines) and the fit with VWA (thick lines). 

• Iterative estimate of T e g, log g, and £: The atmospheric 
parameters and the microturbulence were then refined in 
several steps. This was done by minimizing the correlations 
between the abundance of [Fel/H] lines and equivalent width 
and excitation potential (EP), and requiring good agreement 
between the abundances of [Fel/H] and [Fen/H] (the differ- 
ence in abundance: A(Fei) -A(Feii) is often referred as "ion- 
ization balance"). This task is accomplished by inspecting 
the results on the GUI and calculating the needed models 
as required (see Fig. [5]). 

• Chemical abundance: on the basis of the atmospheric 
parameters determined from [Fel/H] and [Fen/H], the mean 
abundances for all the elements with good lines can be finally 
calculated. 

In Fig. [5] we show the abundances of Fe for the best 
model (T e ff=7700±150 K; log p=4.39 ± 0.06 dex; vsini= 
7.0 ± 1 kms -1 and £=2.4 ±0.2 kms -1 ) and for a model with 
log<? decreased by 0.4 dex, i.e. a value approximately equal 
to that evaluated above with the Mgl triplet. The abun- 
dances are shown versus EW and EP. Open circles are Fel 
and solid circles are Fell lines, respectively. The mean abun- 
dance and rms scatter of Fel and Fell lines are given in each 
panel. We note that in the calculation of the ionization bal- 
ance we adopt the rms of the mean, hence the uncertainties 
on the Fel and Fen abundances are significantly smaller (e.g. 
~ 0.01-0.02 dex). Bearing this in mind, a comparison of top 
and bottom panels in the figure clearly show that in the 
VWA framework a value loggr « 4.0 dex is unlikely. 

To estimate the uncertainties of the atmospheric pa- 
rameters, we repeated the above outlined analysis by vary- 



ing significantly T e g, logg and £ (one parameter is allowed 
to vary while the other two are kept fixed). In this way 
we can determine when the correlations with EW and EP 
become significant or th e ionization balance begins to de- 
viate from equality (see iBruntt. De Cat fc Aertsj 120081 . for 
details). The same perturbed models can be used to es- 
timate the uncertainty on the [Fe/H] value by estimating 
the variation on its value caused by 1 a perturbation of one 
paramet er among T en , logq, £ and tak ing fixed the other 
two (see IBruntt. De Cat fe Aertsl 2008, for a detailed dis- 
cussion). The resulting three uncertainties were summed in 
quadrature together with the rms scatter of Fel to obtain 
the final error on [Fe/H]. We underline that this method for 
the determination of the uncertainties gives only an internal 
estimate since the absolute temperature scale of the model 
atmospheres may be systematically wrong. Thus, our mea- 
sures could show a good precision, but the accuracy is most 
likely not as good (|Bruntt et al.ll2010bl ). 

It is important to note that one of the physical assump- 
tions in the models adopted here is local thermodynamical 
equilibrium (PTE), but deviations from LTE start to become 
important for stars hotter than about 6300 K, especially for 
stars more metal poor than the Sun. Thus we have included 
the NTPE corr e ctions in the present analysis according to 
iRentzsch-Holml (1 19961 ). The correction for neutral iron in 
our case is [Fei/H]NLTE — [Fei/H]LTE + 0.11 dex. When this 
correction is applied, Fen (unaffected by NLTE) must be 
increased by adding +0.2 to logg. This correction has been 
applied in the results from VWA reported here. 

The stellar parameters and the abundances for iron and 
the other measured chemical species obtained by means of 
VWA are shown in Table [3] The uncertainties were calcu- 
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Table 3. Comparison among atmospheric parameters and abun- 
dances derived by spectral synthesis modeling and by VWA ap- 
proach. N denotes the number of lines used with VWA package. 
T e jj is in Kelvin, logg is in dex, while ■usini and £ are in Km/s. 
Abundances are expressed in the form log JV e i/JVn>t ■ An asterisk 
indicate that these uncertainties were rc-dctcrmincd in Sect. 6. 





Synthesis 


V WA 


IN 


Sun 


T c ff 


7500 ± 180 


7700 ± 150 








4.0 ±0.28* 


4.39 ±0.20* 






v sin i 


7.0 ±0.5 


7.0 ±1.0 






£ 


2.4 ±0.5 


2.4±0.2 






C 


-3.81 ±0.12 


-3.85 ±0.08 


7 


-3.60 ±0.05 





-3.29 ±0.10 


— 


- 


-3.34 ±0.05 


Na 


-6.11 ±0.10 


-5.95 ±0.10 


2 


-5.79 ±0.04 


Mg 


-4.62 ±0.14 


-4.77 ±0.15 


1 


-4.43 ±0.04 


Al 


-5.59 ±0.08 


-5.56 ±0.15 


1 


-5.58 ±0.03 


Si 


-4.59 ±0.13 


-4.60 ±0.08 


8 


-4.52 ±0.03 


S 


-5.00 ±0.17 


-4.88 ±0.10 


2 


-4.91 ±0.03 


Ca 


-5.62 ±0.15 


-5.90 ±0.07 


9 


-5.69 ±0.04 


Sc 


_q i q + n i fi 


n nfi + n 06 


4 


_q 88 + 04 


Ti 


-7.14 ±0.10 


-7.20 ±0.07 


21 


-7.08 ±0.05 


V 


-8.28 ±0.21 


-8.33 ±0.15 


1 


-8.10 ±0.08 


Cr 


-6.48 ±0.14 


-6.55 ±0.07 


26 


-6.39 ±0.04 


Mn 


-6.92 ±0.14 


-7.31 ±0.09 


7 


-6.60 ±0.04 


Fe 


-4.68 ±0.15 


-4.72 ±0.06 


125 


-4.53 ±0.04 


Co 


-7.39 ±0.08 


-7.21 ±0.15 


1 


-7.04 ±0.07 


Ni 


-6.13 ±0.08 


-6.02 ±0.06 


26 


-5.81 ±0.04 


Cu 


-8.35 ±0.15 


-8.57±0.15 


1 


-7.84 ±0.04 


Zn 


-7.97 ±0.04 


-7.93 ±0.07 


2 


-7.47 ±0.05 


Y 


-9.86 ±0.11 


-9.82 ±0.09 


4 


-9.82 ±0.05 


Zr 


-9.35 ±0.13 


-9.32 ±0.15 


1 


-9.45 ±0.04 


Ba 


-9.65 ±0.19 


-9.83 ±0.15 


1 


-9.85 ±0.09 


La 


-10.78 ±0.16 






-10.93 ±0.04 


Ce 


-10.46 ±0.19 


-10.70 ±0.15 


1 


-10.45 ±0.04 


Nd 


-10.58 ±0.16 






-10.61 ±0.04 



lated as for iron but only for elements with at least three 
lines measured. For those chemical species for which less 
than three lines have been identified, we arbitrarily set the 
errors in 0.15 dex for elements with only one line and 0.10 
dex for elements with two lines. 



6 COMPARISON OF THE METHODS 

In the previous section we have analyzed the high resolu- 
tion spectrum of the Am: star HD 71297 using two different 
approaches: the classical spectral synthesis and VWA. The 
results for effective temperature, gravity, vsini, microtur- 
bulence, and chemical abundances, obtained with the two 
methods have been summarized in Tab. [3] We can now com- 
pare between them the stellar parameters obtained with the 
two quoted methods, taking also into account the completely 
independent results obtained in Sect. [4] for T c fr and log g as 
summarized in Tab. 

First of all, the atmospheric models we used in our 
calculations are both ATLAS9 based but with different 
convection zone treatments. In particular, in Sect. 15.11 it 
was compute d using the classic al treatment MLT with 
fixed a=1.25 (|Castelli et alJll997h . In Sect. 15. 2| th e mode l 
was interpolated in the fine grid by iHeiter et al] (|2002f ). 
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Figure 7. Upper panel - Comparison between observed and syn- 
thetic H/3 profiles computed with two convection models repre- 
sented by red dashed line (MLT) and by blue dot-dashed line 
(CM). Bottom panel shows the differences between normalized 
fluxes in H/3 profile obtained from two models differing only by 
the convection treatment. 



w hich uses a more adv a nced convection description based 
on lCanuto fe Mazzitellil <|l992h (CM). 

IHeiter et al.l <|20o3 T investigated the effects of the differ- 



ent convection approaches into the Balmer lines profiles and 
they concluded that the observed profiles should have an ac- 
curacy of at least 0.5 % to clearly appreciate the differences 
between convection models with different efficiency. 

In order to quantify these effects on our target, we com- 
puted two synthetic H/3 profiles using the same T e ff = 7700 K 
and log g = 4.39 but with different convection models. We 
show this comparison in Fig. [7] where in the upper panel 
we compare the observed blue wing of the H/3 with overim- 
posed the profiles computed with MLT (red dashed curve) 
and with CM treatment (blue dot-dashed line). In the bot- 
tom panel we report the difference in percentage between 
the two models. Since the maximum difference is of the or- 
der of 1.5%, indistinguishable at our level of signal-to-noise 
ratio, we conclude that for our target we can neglect the 
differences arising from the convection models adopted. 

To verify the accuracy of this conclusion we have re- 
peated the calculation for effective temperature, gravity, and 
abundances as described in Sect. 15.11 Again, we estimated 
T e ff by simultaneous fitting of H7 and H5 but using ATLAS9 
models modified for the CM treatment convection. We ob- 
tained T c ff = 7600 ± 150 K, that is closer to the photomet- 
ric/interferometric value, but also totally consistent with the 
value derived adopting the MLT treatment for convection. 
Using this temperature, we estimated gravity fitting the Mgl 
triplet and we obtain log g = 4.10 ± 0.10, again compatible 
with previous value. 

For what concerns the abundances, the model built with 
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Figure 8. Comparison between abundances computed using 
models that differ for the treatment of the convection theory, 
MLT vs CM. 

these values of effective temperature and gravity has been 
used to repeat the spectral synthesis analysis, adopting the 
same intervals as in Sect. 15.11 In Fig. [8] we displayed the 
differences between the abundances reported in Tab.[3]com- 
puted with MLT based model and those computed with CM 
model. Also in this case the results are consistent, being the 
weighted average of all differences equal to —0.01 ± 0.09. 

Starting with T e g, the classical spectral synthesis 
method and VWA agree well within 1 a. There is also a 
1.3cr agreement with the T cr f value estimated from photom- 
etry and interferometry, albeit the VWA value is closer by 
200 K to the reference value. On the contrary, concerning 
the value of log g, there is a significant discrepancy between 
the classical spectral synthesis method and VWA results: 
log g =4.0 ± 0.1 dex from synthesis of the strong Mgl triplet 
and log g =4.39 ± 0.06 dex from ionization equilibrium of 
Fel/Fen (VWA). The situation is even more complex if we 
take as a reference the gravity evaluated independently from 
spectroscopy in Sect. log <j=4.17±0.05 dex. This value lies 
halfway between the two spectroscopy-based methods, be- 
ing incompatible by more than 1 and 2 a with respect to the 
results of classical spectral synthesis and VWA methods, re- 
spectively. Thus, if we trust in interferometry and parallax 
we have to conclude that even high resolution spectroscopy 
is unable, at least in the case HD 71297, of evaluating the 
logg value better than 0.2 dex, being also the internal errors 
on this value underestimated for both methods. 

To investigate further the accuracy of the results ob- 
tained with the two different spectroscopic analysis about 
T c ff and logg, we built the Spectral Energy Distribution 
(SED) of HD 71297 (see Appendix for all the details on the 
construction of the SED) as shown in Fig. [9] The figure 
shows the comparison of the various photometric or spec- 
trophotometric sources with three models (in different col- 
ors and line styles), as estimated in the previous section^. 
As well known, the Balmer jump is sensitive to surface grav- 
ity (as well as to effective temperature and line blanketing), 
hence looking at the UV part of the SED we can get some 
insight about the "best" couple (T e fF, logg). Even a qual- 
itative check on the figure show that the T e ff, logf/ couple 
evaluated in Sect. 14.11 (from parallax and interferometry) 
reproduces better the SED at all wavelength. In particular, 

7 All the models were calculated for [Fe/H]=— 0.15 as derived for 
HD 71297 



in the far UV, the classical spectral synthesis method has 
the worse agreement due to the low T e s, whereas across 
the Balmer Jump, VWA seems to have a too high logg. 
This reflects the fact that both spectral synthesis and VWA 
methods go through a modeling that introduces errors, then 
it is important to discuss some of them. 

The error we found in logg by fitting the wings of the 
Mgl triplet (i.e. 0.1 dex) is the formal error due to the fit- 
ting procedure. Actually there is at least another source of 
error that we have to discuss and consider. The result of 
the modeling of a spectral line depends essentially on the 
atomic parameters of the transition and in particular, for 
what concern the width of the spectral lines, we have to pay 
attention to the radiative, Stark and Van der Waals damping 
constants. The values of these three numbers set the broad- 
ening of the line and it is straightforward to understand the 
importance of their accuracy on the final value of logg. 

To focus our attention on this problem, we per- 
formed the following simulation: we fitted the observed 
profile of the Mgl A5172.684 A on synthetic profiles by 
varying the damping constants, one at a time, by an 
amount equal to their typical uncer tainties. According to 
iBarklem. Piskunov fc O'Maral (|2000r i, we adopted 5% as 
typical error in damping constants, obtaining that such a 
variation of 7stark and 7waais lead to an uncertainty in log g 
equal to ~ ±0.3 dex. Similarly, a variation of 7 ra( j of 5%, 
lead to a « ±0.15 dex uncertainty on logg. Therefore, con- 
sidering these errors as random and independent, we can 
conclude that the final uncertainties on gravity estimated 
by fitting the wings of one line is given from their quadratic 
sum and then equal to 0.45 dex. However, since we have 
applied this method to all the Mgl triplet, the final uncer- 
tainty coming from atomic parameters is 0.26 dex. Finally, 
considering also the formal error on the fitting procedure, 
we obtained the total error on log g equal to 0.28 dex. With 
such an error, log g derived with classical spectral synthesis 
is consistent with other values found in this study. 

As for VWA, the estimate of log<? by means of the ion- 
ization balance can be affected by over-ionization effects and 
uncertainties in the temperature structure of the model at- 
mos phere, as well as to unc ertainties on the atomic parame- 
ters (|Fuhrmann et alii 19971 ). Thus, the formal error on logg 
we derived in our analysis has to be regarded as a lower 
limit for the true uncertainty on this parameter. Since it 
is difficult to estimate quantitatively the effect of the vari- 
ous source of uncertainty on the final value of loggr, we de- 
cided to assign a "typical er ror" of 0.2 de x . This suggestion 
is justified by the work of iBruntt et~al] (|2012T ) who used 
VWA to analyze high resolution spectra for 93 solar-type 
stars possessing accurate log<? values measured by means 
of asteroseismological techniques on the basis of outstand- 
ing Kepler satellite photometry. These authors find an aver- 
age difference logg(spec) — logg(asteroseismic)=±0.08±0.07 
dex with a few objects showing larger discrepancies (~0.3 
dex). Thus, an uncertainty of the order of 0.2 dex in logg 
is absolutely n ormal. A similar conclusion was reached by 
ISmallevI (120051 ) : "Realistically, the typical errors on the at- 
mospheric parameters of a star, will be T or r±100 K (1~2%) 
for ±0.2 dex (~20%) for logg ". Finally, as for the values of 
iisini, £, and chemical abundances, an inspection of Tab. 
reveals that there is a good agreement between spectral syn- 
thesis and VWA, in spite of the discrepancies in T gf and 
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Figure 9. Comparison between observed spectral energy distribution and four theoretical fluxes. The meaning of the symbols is the 
following: triangles arc TD1 fluxes, filled circles represent Geneva photometry, diamonds arc fluxes from uvby magnitudes, asterisks are 
UBVJHK fluxes, and open circles represent the spectrophotometry. The theoretical distributions have been calculated for: T e jj=7810 K 
and log g= 4.17 (continuous black line), T cff =7700 K and log g= 4.39 (CM model, dotted blue line), T ofi = 7500 K and logg=4.00 (MLT 
model, dashed red line), and T ef f=7600 K and log g= 4. 10 (CM model, dash-dotted purple line) 



logg. The agreem ent in ch e mical abundance is not surpris- 
ing. As argued by ISmallevI (|2005l ). the atmospheric param- 
eters obtained from spectroscopic methods alone may not 
be consistent with the "true" values as obtained by model- 
independent methods, but this is not necessarily important 
for abundance analyses of stars. 



7 CONCLUSION 

In this paper we have tested two different methods com- 
monly used in the recent literature to characterize the 
stellar atmospheres and their chemical pattern: classical 
spectral synthesis and VWA package. This test was per- 
formed analyzing the spectrum of HD 71297 carried out 
with SARG@ TNG. This object has been classified by 
ICowlevI (Il968l) as a suspect ed marginal metallic star, and 
later on Abt fc Levy! (| 19851 ) assigned the spectral type of 
kA8hA9mF0. 

From direct measurements of distance and diame- 
ter, we obtained for our target the following astrophysi- 
cal parameters: T off = 7810 ± 90 K and log g = 4.17 ± 0.05. 
As a by-product we were also able to derive esti- 
mation of R/R = 1.97 ±0.14, M/M = 1.771°; is, and 
age = 790 ± 90 Myrs. 

For what concern the main part of our paper, i.e. the 
comparison of the two different approaches, we can summa- 
rize as follows: 



• Classical spectral synthesis method gives us the follow- 
ing values: T efl = 7500 ± 180 K and log.g = 4.10 ±0.28 that 
are consistent with the previous values, at least within the 
errors. Projected rotational velocity has been evaluated in 



7.0 ±0.5 km s _1 and f = 2.4 ±0.5 km s" 1 . With these pa- 
rameters, abundance analysis gives us a general underabun- 
dance of iron peak elements with the exception of calcium 
that is solar in content. Solar abundances are also shown 
by oxygen, yttrium, zirconium, barium and rare earths. We 
also tested the influence of the convection in the calculation 
of the synth e tic sp ectrum. As in the case of VWA, we used 
iHeiter et al.l (|2002l ) models to repeat the analysis, obtaining 
results totally co nsistent with tho se obtained with MLT the- 
ory as treated in iKuruczl (|l993al ) with no overshooting and 
a = 1.25. Thus we conclude that, at least in the case of a star 
as hot as HD 71297, and for the resolution and SNR of our 
spectrum, the results are only slightly affected by the choice 
about the particular treatment of convection adopted for the 
analys is. The same conclusion is shared by iGardiner et al.l 
(|l999h . who conclude that the MLT and CM models all give 
similarly reasonable results. On the contrary, for cooler tar- 
gets the role of convection will be more significant and we 
will take it into account properly. 

• By using VWA we obtained the following parame- 
ters: Tefl = 7700 ± 150 K and logg = 4.39 ±0.06, 
a projected rotational velocity of 7.0 ±1.0 km s _1 and 
£ = 2.4 ± 0.2 km s _1 . Considering the experimental errors all 
these quantities are comparable with the ones derived with 
the others method. Also the abundances show a pattern sim- 
ilar to the one computed with classical spectral synthesis. 

As a general conclusion we can state that the methods 
considered in this study to derive fundamental parameters 
useful to characterize stellar atmosphere give consistent re- 
sults, if we consider all the sources of experimental errors. 

An important result of this study concern the chemical 
pattern found for HD 71297. Contrary to what is expected 
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from the previous classification, our abundances (reported in 
T ab.[3l) do not look like those of normal Am star (see Fig. 5 
in lCatanzaro et al"1l201ll . for example). In fact, iron-peak el- 
ements show moderate underabundances, as well as carbon 
and sodium. Heavy elements like yttrium, zirconium, barium 
and rare earths, that usually in A-type stars display abun- 
danc es greater than the Sun analogues (|Erspamer & North! 
l2003h . are quite normal in our target. 

The results shown in this paper, concernin the consis- 
tency of classical spectral synthesis and VWA approaches 
to the analysis of stellar spectra, will allow us to confi- 
dently apply these codes to the other Am stars observed 
at SARGQTNG. 
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• JHK magnitudes from 2MASS survey (|Skrutskie et al.l 
20061). converted in physica l un its by using 
van der Bliek. Manfroid fc Bouchetl (|l99rj ). 

As discussed in Sect. 14. ll we neglected the effects of in- 
terstellar reddening. In order to be consistent with the iron 
abundance we found in our previous analysis, each theoret- 
ical SED has to be calculated with an opacity ODF cor- 
responding to a metalicity of [Fe/H] = —0.15. To achieve 
this goal, we computed for each couple of (T g, logg) two 
synthetic fluxes, one for ODF=[0.0] and one for ODF=[-0.5] 
and then we interpolated between them. 



APPENDIX A: SPECTRAL ENERGY 
DISTRIBUTION 

The observed SED has been obtained by merging various 
data sources collected from the literature. In particular, the 
stellar flux shown in Fig. [9] was constructed by using the 
following data: 

• UV fluxes taken from TD1 satellite (| Thompson et all 
1 19781 ) that cover the 1565-2 740 A range; 

• uvby magnitudes from lHauck fc Mermillodl (|l99Sft and 
co nverted in ph ysical units by using the calibrations given 
by iGravl (|l998l ): y magnit ude has been derived fr om UBV 
color via the calibration of iMcSwain fc Giesl (120051) 

• Geneva photometry taken from iRufener (1 19881 ) and 



converted in fluxe s by m eans of the calibrations given by 
IRufener fc Nicoletl (|l988l ): 

UBV magnitudes taken from Mermilliod (Il99ll) and 



converted in fluxes using iBessell. Castelli fc Pled (| 19981) cal- 
ibrations; 

• spectrophotometry in the ra nge 5800-8000 A taken from 
IWeaver fc Torres- Dodgenl (|l995l ); 
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